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Abstract. In this paper, we consider the Universe deep inside of the cell of uniformity. 
At these scales, the Universe is filled with inhomogeneously distributed discrete structures 
(galaxies, groups and clusters of galaxies), which disturb the background Friedmann model. 
We propose mathematical models with conformally flat, hyperbolic and spherical spaces. 
For these models, we obtain the gravitational potential for an arbitrary number of randomly 
distributed inhomogeneities. In the cases of flat and hyperbolic spaces, the potential is finite 
at any point, including spatial infinity, and valid for an arbitrary number of gravitating 
sources. For both of these models, we investigate the motion of test masses (e.g., dwarf 
galaxies) in the vicinity of one of the inhomogeneities. We show that there is a distance 
from the inhomogeneity, at which the cosmological expansion prevails over the gravitational 
attraction and where test masses form the Hubble flow. For our group of galaxies, it happens 
at a few Mpc and the radius of the zero-acceleration sphere is of the order of 1 Mpc, which 
is very close to observations. Outside of this sphere, the dragging effect of the gravitational 
attraction goes very fast to zero. 
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1 Introduction 

According to astronomical observations, there is no clear evidence of spatial homogeneity 
up to sizes ~ 150 Mpc [1] (see also [2, 3]). On the other hand, on the larger scales the 
Universe is homogeneous and isotropic, being well described by the ACDM model [4] which 
is the modern realization of the Friedmann model. It is well known that the Hubble flow, 
i.e. the linear velocity-distance relation for receding motion of galaxies due to the expansion 
of the Universe, is one of the characteristic features of the Friedmann model. The receding 
galaxies form the Hubble flows. Therefore, it seems natural to observe the Hubble flows 
on scales larger than 150 Mpc. On the smaller scales peculiar velocities of inhomogeneities 
(galaxies, group of galaxies, etc) may distort these flows. Using the characteristic values 
of peculiar velocities v ~ 200 -j- 400 km /sec and the present value of the Hubble parameter 
Hq ~ 70 km/sec/Mpc [4], we can get a rough estimate of distances at which the Hubble 
flow velocity is of the same value as peculiar velocities: D = v/Hq ~ 3 4- 6 Mpc. Hence, the 
Hubble flows can be observed at distances greater than this rough estimate. From this point, 
it looks reasonable that Edwin Hubble discovered his law in 1929 after observing galaxies 
on scales less than 20-30 Mpc. However, recent observations [5-9] indicate the presence of 
Hubble flows at distances less than a few Mpc (even less than 3 Mpc (see, e.g., [8])) from the 
center of our group of galaxies. It means that peculiar velocities are sufficiently small to dilute 
considerably the Hubble flow. In other words, the local Hubble flow is cold [8-10]. Obviously, 
these observations require rigorous theoretical substantiation. It is clear that deep inside of 
the cell of uniformity peculiar velocities are not the only reason which can destroy the Hubble 
flow. Here, the gravitational interaction of the inhomogeneities plays an important role, and 
the rough estimate made above does not indicate at which distances from the gravitating 
objects the Hubble flow overcomes this attraction. To resolve this important problem, it was 
suggested in [11-13] that the cosmological constant/dark energy is responsible for the local 
cold flows. This idea was also supported in the papers based on observations [6-8]. However, 
this suggestion was criticized in the article [14], where the authors argue that at present time 
the peculiar velocities are sufficiently strong to dilute the Hubble flows in the local region. 
To prove it, they used the hydrodynamical approach. 

As we have mentioned above, for each inhomogeneity there is a local region where the 
gravitational attraction prevails over the cosmological expansion. Obviously, the character- 
istic size of this region depends on the mass of the inhomogeneity. We can introduce the 
notion of the zero-acceleration surface where a free falling body (which has the zero angular 
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momentum) will have the zero acceleration. Roughly speaking, test bodies undergo the local 
dynamics inside of this surface and follow the global expansion of the Universe outside of it. 
One of the main aims of our paper is to provide a rigorous definition of this surface. To do it, 
we consider our Universe deep inside of the cell of uniformity. At such scales and on late stages 
of evolution, the Universe consists of a set of discrete inhomogeneities (galaxies, groups and 
clusters of galaxies). Hence, classical mechanics of discrete objects provides more adequate 
approach than hydrodynamics with its continuous flows. Obviously, inhomogeneities perturb 
the Friedmann-Robertson- Walker (FRW) metrics. In the weak field limit (the gravitational 
field of the inhomogeneities is weak and their velocities are small with respect to the speed of 
light) such perturbations are reduced to the gravitational potential of the considered system 
(see, e.g., [15, 16]). So, the first problem is to find this potential. In our paper we obtain a 
master equation for the gravitational potential of the system (an arbitrary number of ran- 
domly distributed inhomogeneities) and solve this equation for models with conformally flat, 
hyperbolic and spherical spaces. In the cases of flat and hyperbolic spaces, the potential is 
finite at any point, including spatial infinity, and valid for an arbitrary number of gravitating 
sources. In the case of the flat space, to get such result, we need to cut off the gravitational 
potentials of each inhomogeneity at some concrete distances tq which we call the radius of 
local gravity. The gravitational attraction of the inhomogeneities is absent outside of this 
radius. The cutoff is smooth, i.e. the potential is smoothly matched at tq. Such procedure 
looks a little bit artificial, but the hyperbolic case is free of this procedure. It is worth noting 
that the Schwarzschild-de Sitter solution is valid only in the case of one gravitating mass 
and its potential is divergent at spatial infinity. Another serious drawback of this solution 
from the cosmological point of view is that it does not take into account the average energy 
density of the matter in the Universe, which forms 27% of the total energy density at the 
present stage and depends on time. Our models are free from all of these defects. 



Having the gravitational potential of an arbitrary system of inhomogeneities, we can 
investigate the motion of test masses (e.g., dwarf galaxies) in this field and the formation of 
the Hubble flows by these masses. This approach is very convenient for numerical simulation. 
This is one of the main results of our paper. Additionally, for both of our models, we define 
analytically the radii of the surfaces of the zero acceleration (the zero-acceleration radii). For 
our group of galaxies, they are of the order of 1 Mpc, which is very close to the observations 
[9]. Outside of these surfaces, the dragging effect of the gravitational attraction goes very fast 
to zero. The presence of the cosmological constant is not absolutely crucial for the Hubble 
flows. The reason for the Hubble flow is the global cosmological expansion of the Universe. 
However, the cosmological constant provides the accelerating motion of the Hubble flows. 
Moreover, in the presence of the cosmological constant, the Hubble flow is less diffused by 
peculiar velocities than in its absence (see, e.g., [17]). 

The paper is organized as follows. In section 2 we describe the background model, intro- 
duce the energy-momentum tensor for the inhomogeneities, perturb the background metrics 
by these inhomogeneities and obtain the master equation for the gravitational potential of 
this system. The gravitational potentials as well as equations of motion for test masses form- 
ing the Hubble flows are found in sections 3 and 4 for conformally flat and hyperbolic spaces, 
respectively. The main results are summarized in concluding section 5. 
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2 Scalar perturbations of FRW Universe 



To start with, we consider a homogeneous isotropic Universe described by the Friedmann- 
Robertson-Walker (FRW) metrics 

/ \ 



ds z 



drj 1 — r ) OL pdx a dxP 



drf 



5 a pdx a dx^ 



\ 



1 + -}C (x 2 + y 2 + z 2 ) 



(2.1) 



where tC = — 1,0, +1 for open, fiat and closed Universes, respectively. The Friedmann equa- 
tions for this metrics in the case of the ACDM model read 



and 



A, 



(2.2) 



(2.3) 



where T~L = a'/a = (da/drj)/a and k = 8ttGn/c 4 (c is the speed of light and Gn is the 
Newton's gravitational constant). Hereafter, the Latin indices i,k,= 0, 1,2,3 and the Greek 
indices a, f3 = 1, 2, 3. T lk is the energy-momentum tensor of the averaged pressureless dustlike 
matter 1 . For such matter, the energy density T = pc 2 /a 3 is the only non-zero component, 
p is a constant which we define below. It is worth noting that in the case JC = the 
comoving coordinates x a may have a dimension of length, but then the conformal factor a 
is dimensionless, and vice versa. However, in the cases K, = ±1 the dimension of a is fixed. 
Here, a has a dimension of length and x a are dimensionless. For consistency, we shall follow 
this definition for /C = as well. To make reading easier, we insert the following table where 
explain the dimensions adopted in our paper: 



quantity 


symbol 


dimension 


scale factor 


a 


L (length) 


rest mass density in comoving frame 


P 


M (mass) 


average rest mass density in comoving frame 


P 


M (mass) 



Conformal time r\ and synchronous time t are connected as cdt = adr\. Therefore, eqs. 
(2.2) and (2.3), respectively, take the form 

For the late stages of the Universe evolution, we neglect the contribution of radiation. 
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and 



W-W4 + «a). (2.5) 



a 



6a 3 3 U V 2 a 3 



where ao and -ffo are the values of the conformal factor a and the Hubble "constant" H = 
a/a = (da/dt)/a at the present time t = to (without loss of generality, we can put to = 0), 
and we introduced the standard density parameters: 



— 4 

KpC 



Ac 2 „ Kc 2 



therefore 

n A / + V A + = 1 . (2.7) 
The spatial part of the metrics (2.1) can be also written in spherical coordinates: 

dl 2 = laP dx a dx p = d x 2 + T, 2 ( x )dn 2 2 , (2.8) 

where 

{sinx, X ^ [0, 7r] for K, = +1 

X, XG[0,+oo) for /C = (2.9) 
sinhx, X ^ [0, +oo) for /C = — 1 

In our paper, we consider the stage of the Universe evolution, which is much latter than 
the recombination time. At this stage, the formation of inhomogeneities (stars, galaxies, 
clusters of galaxies) has been generally completed. As we mentioned in Introduction, the 
cell of statistical homogeneity/uniformity size is of the order of 150 Mpc [l]-[3]. On much 
bigger scales the Universe is well described by the ACDM model with matter mainly in the 
form of dark matter plus the cosmological constant. Here, dark matter is well simulated by a 
pressureless perfect fluid, i.e. dust, and the hydrodynamical approach provides the adequate 
description of the model. On scales smaller than the cell of uniformity size, e.g., less than 150 
Mpc, and on late stages of evolution, the Universe is highly inhomogeneous and structured. 
We can see isolated galaxies, which form clusters and superclusters. The observations also 
strongly argue in favor of dark matter concentrated around these structures. Obviously, these 
visible and invisible isolated inhomogeneities can not be represented in the form of liquid. 
Therefore, hydrodynamics is not appropriate to describe their behavior on the considered 
scales. We need to create a mechanical approach where dynamical behavior is defined by 
gravitational potentials. To perform it, we start with the energy-momentum tensor of non- 
interacting (except for the gravity) randomly distributed particles (inhomogeneities, in our 
case) [18]: 

rpik ST^ mpC 2 dx 1 dx k ds sr-^ m p c 2 dx 1 dx k ds 

= (_0)i/2 [ t ] ~dJ~dV^dt d{ - r ~ r p ) = z^ (-5)1/2^] ~dllbd^ d(jc ~ Tp) 

m D c 2 dx 1 dx k dri „ . . 

-4(r - r.) , (2.10) 



E 



v 



(— gY/ 2 [rf\ drj dij ds 



where m p is the mass of p-th inhomogeneity and [t] and [rj] indicate that the determinant is 
calculated from the metric coefficients defined with respect to synchronous t or conformal ij 
times. In the rest of this section we deal with conformal time r/ and we drop for simplicity 
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the symbol [rj]. It is worth noting that for the Universe filled by inhomogeneities, eq. (2.10) 
describes the true energy- momentum tensor of matter at any cosmic scales. 

In the ACDM model, the main contributions come from the cosmological constant and 
the nonrelativistic matter. Therefore, the peculiar velocities should be much less than the 
speed of light: 

dx Q dx a l , <N 

— = a— 2.11 
drj dt c 

Therefore, we can assume that T 00 is the only non-zero component of the energy-momentum 
tensor: 



P 



where 7 is the determinant of the metrics -y a /3 and we introduce the rest mass density 

p = -^^2 m P 6 (. r ~ t p) ■ ( 2 - 13 ) 

After averaging (2.12) over all space, we get T°° = ~pc 2 /a 5 , where p is the average rest mass 
density (2.13) and we use the unperturbed metrics (2.1). Therefore, T = /5c 2 /a 3 as we have 
mentioned above. 

Obviously, the inhomogeneities in the Universe result in scalar perturbations of the 
metrics (2.1). In the conformal Newtonian gauge, such perturbed metrics is [15, 16] 



ds 2 « a 2 



(1 + 2<S>)drj 2 - (1 - 2m) lap dx a dx p ] , (2.14) 



where scalar perturbations $ and * depend on all space-time coordinates 77, x, y, z and satisfy 
equations 

A* - 3W(*' + U$) + = ^a 2 <5T ° , (2.15) 
9 ;{^' + U<S>) = Iko, 2 6T2 = 0, (2.16) 



dx? y ' 2 



-- K a 2 5T$ = 0, 



(2.17) 



where the Laplace operator A and the covariant derivatives are defined with respect to the 
metrics 7^ . The condition 5T® = follows from the nonrelativistic nature of the considered 
matter, i.e. \8Tp\ <C 5Tq, and we can drop 5T® with respect to 5Tq. To clarify this point, 
we want to stress that according to Eqs. (2.15) and (2.16), both 5T$ and 5To contribute 
to the gravitational potential $. However, due to the condition (2.11), peculiar velocities 
of inhomogeneities are nonrelativistic and the contribution of 5T® is negligible compared to 
that of 5Tq. In other words, account of 5T® is beyond the accuracy of the model. This 
approach is fully consistent with [18] where it is shown that the nonrelativistic gravitational 
potential is defined by the positions of the inhomogeneities but not by their velocities (see 
Eq. (106.11) in this book). In the case of an arbitrary number of dimensions, a similar result 
was obtained in [19]. The perturbed matter remains nonrelativistic (pressureless) that results 
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in the condition STg = 0. Below, we define the accuracy, with which the latter statement 
holds. Following the standard argumentation (see, e.g., [15, 16]), we can put <I> = ^, then 
the system of above equations reads 

A$ - 3W($' + n$) + 3/C$ = ^a 2 5T ° , (2.18) 
9 r (S' + U§) = 0, (2.19) 



From eq. (2.19) we get 



$" + 3ft $' + (2ft' + ft 2 ) $ - K$ = . (2.20) 

<S>( v ,r) = ^- (2.21) 

where ip(r) is a function of all spatial coordinates and we have introduced c 2 in the denom- 
inator for convenience (for such normalization, if has the dimension (c 2 ) x (length)). The 
function $ is defined up to an arbitrary additive function f(rj) which does not depend on the 
spatial coordinates. Therefore, this function / is not related to the inhomogeneities and we 
can drop it. Below, we shall see that </?(r) ~ 1/r in the vicinity of an inhomogeneity, and the 
nonrelativistic gravitational potential <&(r], r) ~ l/(ar) = 1/R, where R = ar is the physical 
distance. Hence, $ has the correct Newtonian limit near the inhomogeneities. Substituting 
the expression (2.21) into eq. (2.18), we get the following equation for ip: 

A(p + 3/Cv? = ^Kc 2 a 3 5T . (2.22) 

The left hand side of this equation is independent of 77, hence, the right hand side also should 
not depend on rj. This is possible only if 5T$ ~ l/o 3 . Let us now define the conditions of 
implementation of this statement. According to eq. (2.12), we have 

T o = VT^p {2 23) 

Taking into account the perturbed metrics (2.14), we get in the first approximation 



n o n o n o n <± 



where 5p is the difference between real and average rest mass densities: 

5 p = p-p. (2.25) 

In the right hand side of eq. (2.24), the second term is proportional to 1/a 4 and should be 
dropped because we consider the nonrelativistic matter. This is the accuracy of our approach. 
Hence, the perturbation of the energy-density reads 

«5T ° = ^. (2.26) 

Finally, eq. (2.22) takes the form 

Atp + 3}Cip = 4irG N (p-p) . (2.27) 
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This is our master equation where the appearance of p plays a crucial role. It is worth 
mentioning that ip is not a "physical" gravitational potential, and eq. (2.27) looks only 
formally as the usual Poisson equation. Here, p and ~p are comoving local and average rest 
mass densities, respectively, which do not depend on time. The connection between <p and 
the physical potential $ is given by eq. (2.21). 

We now turn to eq. (2.20). Taking into account the relations 

= _J£L = _ft$, $" = -n'$ + n 2 $, (2.28) 

we can rewrite this equation in the following form: 

(W -H 2 - 1C)$ = 0. (2.29) 

On the other hand, eqs. (2.2) and (2.3) show that 

2(W-U 2 -K) -o 1 , . 

-i - 2 = -kT° ~ - . 2.30 

Therefore, for $ from (2.21), eq. (2.20) is satisfied with the adopted accuracy O (l/a 4 ). 

It is worth noting that for an arbitrary perfect fluid with the equation of state p = cue, 
eq. (2.30) reads 

2 {%' - H 2 - K) 



,2 



-K(e + p) = -K(l + w)e. (2.31) 

a* 

For the cosmological constant, the right hand side is equal to zero. However, in the case 
uj 7^ —1, the right hand side behaves as l/a 3 ^ 1+w - ) . The demand e ~ 0(1/ a n ), n > 3 leads 
to the condition co > 0. Therefore, homogeneous quintessence (—1 < oj < 0) and phantom 
matter (to < —1) can not be an alternative to the cosmological constant. In the papers 
[20, 21], it was also pointed out that the quintessence has to be inhomogeneous. 

Coming back to the perturbed metrics (2.14), we can write it now in the form 

ds 2 » (1 + 2$) c 2 dt 2 - a 2 (1 - 2$) laP dx a dxP . (2.32) 

We shall use this metrics for investigation of motion of nonrelativistic test masses. The action 
for a test body of the mass m can be written in the following form [18]: 

S = -mc J ds^-mc J {(l + 2$)c 2 - a 2 (1 - 2<S>)v 2 } l/2 dt , v 2 =^x a x p , (2.33) 

where v is a comoving peculiar velocity which has the dimension (time) -1 . Hence, the 
corresponding Lagrange function has the form 

,2 J 1 „2« 2 / 1 o,m1 V2 ~ A , <P a 2 v 2 \ _ 2 mp ma 2 v 2 



L « -mc 2 n + 2* - a 2 ^r(l -2$)\ » -mc 2 1 + ^- - — — = -mc 2 ^+ 

[ c z J \ ac z 2c z J a 2 

(2.34) 

where we dropped the term O (l/c 2 ). For nonrelativistic test masses, we can also drop the 
term mc 2 . To get the Lagrange equations, we need now the form of the gravitational potential 
ip which is the solution of eq. (2.27). Below, we investigate eq. (2.27) separately for the cases 
K = and K, = ±1. 



-7- 



3 Hubble flows in conformally flat space (/C = 0) 



In this case eq. (2.27) reads 

A<p = 4xG N (p-p), (3.1) 

where the Laplace operator A = Yla=i & 2 1 (dx a ) 2 and the rest mass density p is defined 
by (2.13) where 7 = 1. Obviously, the presence of p destroys the superposition principle 
(it is hardly possible to single out the contribution of each of the inhomogeneities to the 
average density ~p). There is also the well known problem related to the solution of the 
Poisson equation for the infinite space with homogeneous energy density (the Neumann- 
Seeliger paradox) or infinite number of gravitating masses. To find a finite solution of such 
equation, we need to assume the form of the spatial distribution of these masses. Therefore, 
to avoid this problem in our case, we suppose that in the vicinity of each inhomogeneity, the 
gravitational potential is defined by the mass of this inhomogeneity and is not affected by 
other masses. We consider a simplified version where the inhomogeneities are approximated 
by point-like masses, which do not interact gravitationally with each other. We can easily 
generalize this picture to the case where some of inhomogeneities form a gravitationally 
bound system. In this case, we consider such system as one point-like mass concentrated in 
the center of mass of this system. Further, we assume that each point-like mass moi (here, we 
introduce the subscript to differ these gravitating masses from a test mass m) is surrounded 
by an empty sphere of the radius roi defined below and this sphere, in turn, is embedded in a 
medium 2 with the rest mass density ~p (see Figure 1). Below, we shall demonstrate that such 
supposition of the spatial distribution of matter provides the finiteness of the gravitational 
potential at any point of space and for an arbitrary number of inhomogeneities. 
Therefore, inside of each of these spheres, the Poisson equation (3.1) reads 

1 d 2 

£<p = -- r ^(r^)=4irG N [m 8(r)-p], r < r , (3.2) 
r ar z 

where, for simplicity, we omit the subscript % for the mass mo, the radius ro and the gravi- 
tational potential ip. In the cosmological medium, the Poisson equation is 

Av? = 0, r>r . (3.3) 

For this equation we choose the trivial solution 

tp = ^ = , r > r . (3.4) 

ar 

It can be easily seen that the solution of eq. (3.2) with the Newtonian limit at r = 0, which 
is smoothly matched with (3.4) at r = ro, is 

G N m G N m 2 3G N m ^ . . 

r 2tq 2ro 



2 Obviously, in the region of space where the energy density coincides with the average energy density 
pc 2 /a 3 , we restore the unperturbed Friedmann Universe with the ideal cosmological medium characterized 
by p. This ideal cosmological medium takes into account the average effect of all inhomogeneities. It seems 
tempting to solve the problem of missing dark matter in the local Universe with the help of this cosmological 
medium which may play the role of a "smooth ocean" of dark matter suggested in [9]. As we shall see below, 
test masses outside of the spheres (i.e. for r > ro) form Hubble flows which are not affected by the gravitational 
attraction of the inhomogeneities. The gravitational field of an inhomogeneity stops to act outside of this 
sphere. For this reason, we refer to these surfaces as the spheres of local gravity, and the radius ro is the 
radius of local gravity. 
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Figure 1. For the model K, = 0, we depict schematically an example of the large-scale structure 
with five gravitationally unbound inhomogeneitics (e.g., groups of galaxies). Each of inhomogencities 
is surrounded by the vacuum sphere of local gravity with the Poisson equation (3.2) inside of this 
sphere. The radius roi of the sphere and the mass moi of the inhomogeneity are related by eq. 
(3.6). All these spheres are embedded into ideal cosmological medium (the shaded region), where the 
gravitational potential is absent. Inside of this region, test masses follow the Hubble flows diluted by 
peculiar velocities (see eq. (3.13)). For definiteness, we consider the spheres of local gravity of radii 
?'oi = r ,r 2 = 2r , r 03 = 3r , r oi = 4r and r 05 = 5r . 



where the matching condition gives 

Eq. (3.6) shows that ro is just the radius of a sphere of the rest mass density ~p and the total 
mass mo = 4-7rrQp/3. The second term in the right hand side of eq. (3.5) is divergent at large 
distances: 

27rG N p 2 G N m 2 

r = g — r — > +oo lor r — > +oo . (3.7) 

3 2r 

This was the reason of the cutoff of the potential (3.5) at r = tq. The nonrelativistic 
gravitational potential in the Schwarzschild-de Sitter model is not free from the similar 
quadratic divergence. It is worth mentioning that ro and 7j are comoving quantities which 
do not depend on time, ro is dimensionless, ~p has the dimension of mass and <p has the di- 
mension (c 2 ) x (length). The physical rest mass density is ~p ph = p/a 3 . It makes sense 
to estimate at present time (a = ao) the physical radius of the sphere of local gravity 
Ro(a = ao) = Rq = ao^o = ao [3mo/(47rp)] 1 ^ 3 = ^Gnttiq/ (HqQj\j)^ for our group of 
galaxies with mo ~ 1.9 x 10 12 Mq 3.8 x 10 45 g [9]. According to the seven-year WMAP ob- 
servations [4], H w 70km/sec/Mpc « 2.3 x 10~ 18 sec _1 , M ~ 0.27 and A ~ 0.73. Hence, 
^o ~ 7.1 x 10 24 cm » 2.3 Mpc. 

Now let us investigate motion of a test mass m in the considered background. The equa- 
tion of motion is defined by the Lagrange function (2.34) where the gravitational potential 
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is given by eqs. (3.4) and (3.5). The Lagrange equation is 



at a or 

where for the flat comoving space v = dr/dt. According to eq. (3.4), in the cosmological 
medium, the gravitational potential ip = and integration of (3.8) gives 

t 

ci f 1 

v = -z, r = ci / —zdt + c 2 , (3.9) 
a z J a z 

to 

where ci and c 2 are constants of integration defined with respect to the initial time tQ. v is the 
peculiar velocity in the comoving coordinates. The physical distance is R = ar. Therefore, 
the physical velocity V = dR/dt = d(ar)/dt is related to v as follows: 

V = -R + av, (3.10) 
a 

This equation shows that the physical peculiar velocity is v p h = av. With the help of the 

relation (3.10), we can rewrite the Lagrange function (2.34) with respect to the physical 
quantities V and R: 

9 mtp mV 2 md 2 R 2 md 

L= - rac 2-_^ + _ + _-_ ( VR). (3.11) 

The corresponding Lagrange equation is 



— V--R = — -^ + ^R--V. 3.12) 
dt V a I a an a z a 



In the case (p = 0, i.e. in the cosmological medium (r > ro), we can easily integrate this 
equation: 

t 

V = -R+— , R = oci / 4rdt + ac 2 . (3.13) 
a a J a 

to 

Therefore, in the absence of the gravitational potential, the physical velocity consists of two 
parts. They are the Hubble velocity which is "diluted" by the peculiar velocity. Obviously, 
for if = 0, the Lagrange equation defines the Hubble flows. Here, eq. (3.12) in the case of 
radial motion can be rewritten in the form 3 

where we use the Friedmann equation (2.5). It can be easily seen that in the cosmological 
medium (r > ro) the Hubble flows exist even in the absence of the cosmological constant. It 
is therefore not correct to say that the Hubble flows at small distances (a few megaparsecs) 
are due solely to the cosmological constant (see, e.g., [13]). The reason for the Hubble flow is 
the global cosmological expansion of the Universe. However, the acceleration {dV/dt > 0) is 



3 



For radial motion, the equation aR = dR =>■ d (a 2 r) /dt = is equivalent to eq. (3.8) if tp = 0. 
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possible only in the presence of the positive cosmological constant. The cosmological constant 
also reduces smearing of the Hubble flow by peculiar velocities [17]. 

Obviously, in spherical coordinates r,0 = Tr/2,ip the Lagrange function (2.34) reads 

L = -mc 2 -^ + ^fr 2 +rV N ) • (3.15) 
a 2 V / 

The corresponding Lagrange equations are 

^ (ma 2 r 2 ip) =0 => ip = — ^— (3.16) 

at V / ma z r z 

and 

d , a 2+\ _ 1 gy , m 2 (317) 

dt ' a dr m 2 a 2 r 3 ' 
where M = const is the angular momentum. Let us investigate the latter equation in the 
vicinity of the gravitating mass ttiq, i.e. in the region r < where the gravitational potential 
is given by (3.5): 

± ( a 2 r) = --( GNm ° - GNm \ ) + M2 . (3.18) 
dt a \ r 2 r o / m 2 a 2 r 3 

After some algebra where we use the Priedmann equation (2.5), this equation takes the form 

It is interesting to note that this equation contains the cosmological constant A, rather than 
the average density p in contrast to the expression for the gravitational potential (3.5). Eq. 
(3.19) can be also obtained from the Lagrange function (3.11) with (p from (3.5). Taking into 
account the Friedmann equations (2.4) and (2.5), this function reads 

Givmnm m Ac 2 „ 9 mV 2 9 - 9 „ 9 - 9 „ „ 
L= N R ° +y— R + — ■ V 2 = R 2 + R 2 ^j 2 , R<R . (3.20) 

In classical mechanics, this Lagrange function describes motion of a test mass m in the 
gravitational field with the potential 

_ G N m 1 Ac 2 2 d(p _ G N m Ac 2 
V- ^"2—^' -8R-— + — R > R ~ R °- (3 ' 21) 

Exactly this potential enters in the Schwarzschild-de Sitter metrics: 

ds 2 =( 1 + ^j) c2 ^ 2 - f 1 + ) <*fl 2 - ^ 2 (^ 2 + sin2 dd ^ 2 ) ■ ( 3 - 22 ) 

As it follows from (3.21), the derivative —dCpjdR equals at R = R* = [3GNmo/ (Ac 2 )] 
Therefore, Rq/R* = (2n A /n M ) 1/3 (a/a ). Since at present time (a = ao) the parameters 
SI A and fiju satisfy the inequality J7a ~ 0.73 > CIm « 0.27, we get that i?* < /?o (where 
i?o — -Ro( a — a o))- Hence, inside of the sphere of local gravity, the attraction dominates 
for R < R* and repulsion dominates for R* < R < Rq. If M = 0, then R =0 (see eq. 

(3.19)), and for the model AC = 0, R* can be associated with the radius of the zero-acceleration 
surface. In some articles this surface is called the zero-velocity surface [8, 9, 22, 23] or the 
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surface of zero gravity [11]. According to the observations [8, 9, 23], the radius of this surface 
is 0.96 Mpc for our group of galaxies. The calculated value R* = Rq (20a/^a/) -1 ^ 3 ~ 1-3 
Mpc is very close to the observed one. 

Taking into account that Gj^tuq/Rq = A-kGnJ)/ (3a 3 ), it can be easily seen that at 
R = Ro eq. (3.19) is matched to eq. (3.14) (for radial motion M = 0). Outside of the sphere 
of local gravity (R > Ro), test masses form the Hubble flows (3.13). Because $7 a > £2m, 
these flows are accelerated at the late stage of the Universe evolution (see eq. (3.14)). 

4 Hubble flows in conformally hyperbolic space (/C = —1) 

At the beginning of this section, we consider the general case of nonzero curvature (K. = ±1) 
and then dwell on the physically more interesting case of the hyperbolic space (IC = — 1). 
To start with, we turn back to our master equation (2.27). For IC ^ 0, we can define a new 
function (f>: 

4ttG n p 

4> = v+ 3jC , (4-1) 
which satisfies the Helmholtz equation 

+ 3IC(j) = AitGn p , (4.2) 

where p is defined by eq. (2.13). It can be easily seen that the principle of superposition 
holds for eq. (4.2). Therefore, we can solve eq. (4.2) for one gravitating mass. Then, the 
solution for an arbitrary number of masses is the sum of the solutions for each separate mass 
(see below). So, let us consider one gravitating mass which we denote as mo- Without loss of 
generality, we may put this mass in the origin of coordinates. Because of spherical symmetry 
of the problem and for the metrics (2.8), eq. (4.2) in vacuum reads 



1 d 



£%)-^ +3/C0 = O, (4.3) 



£ 2 (x) dx V d x 
where X(x) is defined in (2.9). 

First, we consider the case of the spherical space IC = +1 when X(x) = sinx- The 
solution of this equation with the Newtonian limit at x = is 



4> = 2Ci cos x ~ G N m 2 sin x) ■ (4.4) 

\smx ) 

For any value of the constant of integration C\, this solution (as well as its derivative —d(j)/dx) 
is divergent at x = tt- This reflects the fact that the surface area of the sphere which surrounds 
the gravitating source shrinks to zero for % — > ir. The procedure of the cutoff which we explore 
in the flat space does not work here because equations <p(xo) = an d d(p/dx(x = Xo) = 
are incompatible for Vxo £ [0,7r]. Therefore, we disregard this model. 

Now, we consider the hyperbolic model with /C = —1 where £(%) = sinhx. It can be 
easily seen that in this case the solution of eq. (4.3), which has the Newtonian limit at x ~ > 0: 
<f> — > — Gatttt-o/x an d is finite at infinity: <f> — > for x - > +oo, reads 

, n exp(-2x) . , 

4> = -G N m — — . (4.5) 

smnx 
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This formula demonstrates a number of advantages with respect to the fiat space case. First, 
the presence of the exponential function enables us to avoid the gravitational paradox. In 
some models, such exponential function was introduced by hand. In our model, it appears 
quite naturally. So, we do not need to introduce the sphere of cutoff of the gravitational 
potential. The gravitating masses in this case can be distributed completely arbitrarily. 
Second, we can apply the principle of superposition. For example, the function (f> for all 
gravitating masses is 

exp(-2/j) 



-G N y^mQj- 



sinh L 



(4.6) 



where Zj denotes the geodesic distance between the i-th mass moi and the point of observation. 
Then, according to eq. (4.1), the gravitational potential of this system is 



^ exp(-2Zi) AnG N p 



(4.7) 



To investigate the Hubble flows in this model, we consider motion of a test mass m 
in the vicinity of the gravitating mass uiq where we can neglect the effect of other distant 
gravitating inhomogeneities, i.e. i = 1 in the formula (4.7). In other words, the function <f> 
is given by eq. (4.5). It is easy to verify that Lagrange equations in the equatorial plane 
9 = ir/2 for the considered case have the form 



— (ma 2 sinh 2 yijj 
dt V AY 







1> 



M 



ma 2 sinh x 



and 



cZ / 9 . ■> 1 ocp 
-7; (a \) =— 7T + 
cZi a ox 



M 2 cosh x 
2 a 2 sinh 3 x 



m 



(4.8) 



(4.9) 



It makes sense to estimate x f° r the scales relevant to our model, i.e. for R < 150 Mpc. 
In the hyperbolic space, the physical distance is R = ax- Obviously, R <C a => \' « 1. 
More precise limitation can be found if we know the value of the scale factor ao at present 
time which we can get from the density parameter 0^ = _i = c 2 / (gsq-Z/q) . According to the 
seven-year WMAP observations, Qjc=-i < 8.4 x 10 -3 (see section 4.3 in [4]). Let us take 
VL K= _i ~ 1CT 4 . Then, for H w 2.3 x 10 _18 sec _1 , we get a «4x 10 5 Mpc. Therefore, x < 1 
for all distances R <C 10 5 Mpc, and here we may approximate the function <\> as follows: 
(f) ~ —Gntuq/x- For such values of x, the Lagrange equation (4.9) can be rewritten in the 
form 



R= a -R 
a 



G N m M 2 



R 2 



m 2 R? 



±n M ^ + n A )R 



G N m 
R 2 



+ 



M 2 



m 2 R? 



(4.10) 



where we used the Friedmann equation (2.5). It is clear that the first term on the right 
hand side of this equation is responsible for the Hubble flow (see, e.g., eq. (3.14)). However, 
the second term (the gravitational attraction) works in the opposite direction trying to bind 
gravitationally the system. It makes sense to estimate the distance where both of these terms 
are of the same order: 



R 



G N m 
R 2 



R' 



1 a=a,Q 



R 



G N m 



H 



\d/a\ 



G N m 



a=ao 



H 2 \n A - Q M /2\ 



(4.11) 
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Therefore, at present time, the Hubble flows begin to prevail over the gravitational attraction 
starting from Rh- Obviously, Rh plays the role of the radius of the zero-acceleration surface. 
It can be easily seen that the radius of the local gravity sphere in the flat space case Rq = 
[2G N m /(Hfo M )] 1/3 is related to R H as follows: R H = Rp X |(20 a _/^m) - 1|" 1/3 - For 
f&A ~ 0.73 and fijf ~ 0.27, we get for our group of galaxies Rh ~ 0.61 Rq ~ 1.4 Mpc, which 
is rather close to the observed value 0.96 Mpc [9]. 

For R > Rh, the Hubble flows are distorted by the gravitational attraction (it is the 
so called dragging effect [8, 9]). To quantify these distortions, it is convenient to rewrite eq. 
(4.10) in dimensionless units: 

n R M 2 

g=-, t = H t, M 2 = - _ 4 ■ (4.12) 

°o Rh H$m 2 R H 



Then, we get 



d2R -( n " +nA W" A -> /2 + g. (4.13) 



dP V 25 3 A J R? R3 



Without loss of generality, we may consider the case of the radial motion: M = 0. Obviously, 
at the stages of the accelerating expansion of the Universe, it is reasonable to study the Hubble 
flows at distances, where the right hand side of this equation is non-negative: 

i.e. R > 1 =/- R > Rh at the present time (a = 1). Then, we can easily calculate the 
distance, at which the contribution of the absolute magnitude of the gravitational attraction 
to the right hand side of eq. (4.13) is n% . 



R 



71% 



ioo% \ n A -n M /2 1 1 / 3 
n% J n A -n M /(2^)\ ■ { ' 



For example, at present time, R\q% ~ 2.2, R 5 % 2.8 and R\% ~ 4.7. For our group 
of galaxies it gives, respectively, R\o% ~ 3.1 Mpc, R 5 % 3.9 Mpc and R\% w 6.6 Mpc. 
Therefore, in this model, the Hubble flows are observed at sufficiently small distances much 
smaller than the cell of homogeneity size. To confirm this conclusion, we consider now some 
plots of the velocity V = dR/dt of a test body. In dimensionless units it reads 

y = dR = = ]_y (416) 

dt R H H 

For the Hubble flow V = H(t)R = H{t)R H R we get 

V = ^Mr. (4.17) 
Ho 



4 At the late stage of the Universe evolution, the right hand side of eq. (4.13) (where for simplicity we 
put M = 0) consists of positive and negative terms and is equal to zero at present time (a — 1) for R = 1. 
In this case n% runs from to infinity: < n% < +00, which can be easily seen from the expression 

n% = (r 3 - l) 100% , R > 1 obtained for 5 = 1. 
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Hence, V = R only at the present time t = to. Without loss of generality, we usually put 
t = => t = 0. 

On Figure 2, we depict the ratio V/R for two cases (in both of them M = 0). The 
upper (red) line is a result of numerical integration 5 of eq. (4.13) for the pure Hubble flow 
case, i.e. when the second term (the gravitational attraction) in the right hand side of this 
equation is absent. The middle (blue) line demonstrates how the gravitational attraction 
changes the Hubble flow, i.e. this is the ratio V/R in the presence of the gravitational 
attraction. The initial (at t = to = 0) conditions are taken as follows: R = 1 =4> R = Rh 
and V/R = 1 => V/R = Hq. The dashed line is the eleven-percent deviation from the 
Hubble flow. Our calculations show that the blue line is entirely within this deviation. 



1.00 



0.95 J 



0.90 - 



0.80 




Figure 2. Here, we depict the ratio V/R from the numerical solution of cq. (4.13) without 
gravitational attraction (the upper/red line), i.e. in the case of the pure Hubble flow, and in the 
presence of this attraction (the middle/blue line). The value t = corresponds to the present time. 
The effect of gravitational attraction does not exceed the eleven-percent deviation (the dashed line) 
from the Hubble flow. 



On Figure 3, we depict the evolution of the velocity V of a test mass with the distance 
R in the case of the zero angular momentum M = 0. Solid and dashed lines correspond to 
presence and absence of the gravitational attraction, respectively. All lines start at R = 1 =>■ 
R = Rh with the following initial velocities: V = 1 (blue lines), V = 2 (green lines) and 
V = 3 (red lines). These graphs show that, first, the effect of the gravitational attraction is 
not strong for R > Rh and, second, all lines go asymptotically to the dashed blue line which 
describes the pure Hubble flow (i.e. the admixture of peculiar velocities is absent and the 
test body moves solely due to the expansion of the Universe). Peculiar velocities result in 
deviation from the pure Hubble flow, and the plots demonstrate how fast test bodies with 
different initial values of V join this flow. For our group of galaxies, there is observational 



5 To perform such integration, we need to know the time dependence of the scale factor a(t) which satisfies 
the Friedmann equation (2.4). Obviously, because 0.x. <C 1, we can drop the curvature term, and the solution 
of this equation is common for any value of KL. Then, in dimensionless units, this equation reads (a~ 1 da/dt) 2 = 
Flui/a 3 + ^a- The solution of this equation (with the boundary condition 5 = 1 at t = to = 0) is 



1/3 



1 + 



1/2 



sinh I ^-£l 1 J 2 i 



1/2 



cosh ( ^Q. 1 / 2 ! 



2/3 
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evidence [8, 23] that the peculiar velocities are of the order of 30 km/sec within 3 Mpc, i.e. 
the initial value V ~ 1.3. It can be easily shown that in this case the corresponding solid 
line coincides with the pure Hubble flow already at i? = 2 => fi«2,8 Mpc. 




Figure 3. This plot shows the evolution of the velocity V of a test mass with the distance R. All 
lines start at R = 1 =>• R = Rh with the following initial velocities: V = 1 (blue lines), V = 2 
(green lines) and V = 3 (red lines). Solid and dashed lines correspond to presence and absence of 
the gravitational attraction, respectively. The blue dashed line describes the pure Hubble flow. For 
our group of galaxies, dimensional velocities and distances are V « V x 100 km/sec and Rm Rx 1.4 
Mpc. 

On Figure 4, we depict the effect of the angular momentum M on the motion of test 
masses. Solid and dashed lines correspond to presence and absence of the gravitational 
attraction, respectively. All colored lines start at R = 1 =^ R = Rh with the zero initial 
velocity V = and have the angular momentum M = 1 (blue lines), M = 2 (green lines) and 
M = 3 (red lines). Black solid and dashed lines start with the initial velocity V = 1 at R = 1 
and have the zero angular momentum, i.e. the dashed black line here describes the pure 
Hubble flow. First, these colored graphs demonstrate that test masses with the zero radial 
component of the velocity escape from the zero-acceleration surface outward due to the effect 
of the centrifugal force. Second, the radial components of the velocity go asymptotically to 
the values for the pure Hubble flow. The smaller the angular momentum is, the sooner it 
happens. For example, the blue lines with M = 1 coincide with the corresponding black lines 
at R sa 2. 

It is also of interest to draw the characteristic distances Rq = ^Gnitiq/ (HqQ,m)] , 
-R* = Rq {2VLx/^m)~ 1 and Rh = Ro x \{^a/^m)— 1|~ 1//3 as the functions of the gravitating 
mass rriQ. Therefore, we depict these graphs on Figure 5. Here, solid black (upper), blue 
(middle) and red (lower) lines correspond to Rq, Rh and i?*, respectively. These distances 
are given in Mpc. Vertical dashed orange (left) and green (right) lines correspond to our Local 
Group of galaxies (m « 1.9 x 10 12 Mq [9]) and the Virgo cluster (m « 1.2 x 10 15 Mq [24]), 
respectively. This plot shows that at the present moment the radii of the zero acceleration 
sphere Rh for /C = — 1 model and R* for K. = model are very close to each other. 

To conclude this section, we should note the following. It is obvious that the contri- 
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Figure 4. This plot shows the effect of the angular momentum M on motion of test masses. Solid 
and dashed lines correspond to presence and absence of the gravitational attraction, respectively. All 
colored lines start at R = 1 => R = Rh with the zero initial velocity V = and have the angular 
momentum M = 1 (blue lines), M = 2 (green lines) and M = 3 (red lines). Here, V is the radial 
component of the velocity of a test body. Black solid and dashed lines start with the initial velocity 
V = 1 at R = 1 and have the zero angular momentum, i.e. the dashed black line describes the pure 
Hubble flow. 



bution of the constant term in the gravitational potential (4.7) does not affect the motion 
of nonrelativistic objects (see, e.g., eq. (4.9)). However, this term contributes also in the 
metrics (2.32) and at present time is of the order of 8ttGn7)/ (3aoc 2 ) = QmHqO^/c 2 ~ 2x 10 3 
where ao ~ 4 x 10 5 Mpc. This large value does not contradict our approach because we need 
to take into account the gravitational contribution (4.6) of all inhomogeneities. To see it, let 
us average the i-th component in the first term of (4.7) over some finite volume V: 



— 47T 



XO 



-G N m . 



exp(-2x) 
sinhx 



sinh xdx 



4irG N m i 
3V 



Then, letting the volume go to infinity (xo 
all gravitating masses, we obtain 



+oo 



V 



'total 



1-2 

2 



-xo 



_ e -3X0 



(4.18) 

+oo) and taking into account 



(4.19) 



where p = lim ^m-oi/V. Therefore, the averaged gravitational potential (4.7) is equal to 

zero: Tp = 0. By the way, this is one more advantage of the considered model with tC = — 1 
as compared to the model with 1C = where it can be easily shown that the gravitational 
potential averaged over all spheres of local gravity has a non-vanishing negative value. 
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Figure 5. This plot shows the characteristic distances Rq (black upper line), Rh (blue middle line) 
and i?* (red lower line) as the functions of the gravitating mass tuq. Vertical dashed orange (left) and 
green (right) lines correspond to our Local Group of galaxies and the Virgo cluster, respectively. 

5 Conclusion 

In this paper, we have considered our Universe at scales much less than the cell of homo- 
geneity size which is approximately 150 Mpc. At such distances, our Universe is highly 
inhomogeneous and averaged Friedmann approach does not work here. We need to take 
into account the inhomogeneities in the form of galaxies, groups of galaxies and clusters of 
galaxies. All of them perturb the FRW metrics. We have investigated these perturbations in 
the weak-field limit where the 1/c 2 correction term in the metric coefficient goo defines the 
gravitational potential. The main goal was to get the expressions for this potential and to 
investigate the dynamical behavior of test bodies (i.e. the dwarf galaxies) in the field of this 
potential. First, we have found the master equation (2.27) for the potential <p(r) which is 
conformally related to the physical gravitational potential §(r],r). Formally, this equation 
has the form of the Poisson equation. A distinctive feature of this equation is the presence 
(in the right hand side) of the average comoving rest mass density of matter in the Universe. 

Then, to solve this equation, we have studied in detail two types of models, namely with 
conformally flat K, = and conformally hyperbolic /C = — 1 spaces. These models have both 
similarities and significant distinguishing features. For both of our models, gravitational po- 
tentials are finite at any point, including spatial infinity 6 , and these solutions are valid for an 
arbitrary number of gravitating masses. Note that the Schwarzschild-de Sitter solution was 
found for one gravitating source and the corresponding gravitational potential is divergent at 
infinity. For both of these models, we have shown that at present time there is a distance from 
inhomogeneities at which the cosmological expansion prevails over the gravitational attrac- 

6 In the case of the spherical space /C = +1, the gravitational potential is divergent at \ = 7r - 
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tion. For our group of galaxies the radii of the zero-acceleration sphere (R* for fC = model 
and Rh for fC = — 1 model) are of the order of 1 Mpc which is very close to the observations 
[9]. Outside of these spheres, the dragging effect of the gravitational attraction goes very fast 
to zero. We have also shown that the presence of the cosmological constant is not absolutely 
crucial for the Hubble flows. The Hubble flows may also take place in the absence of the 
cosmological constant. The reason for the Hubble flow is the global cosmological expansion 
of the Universe. However, the cosmological constant provides the accelerating motion of the 
Hubble flows. Moreover, it makes sense to talk about the zero-acceleration surface only if at 
present time > £Im/2- Additionally, the cosmological constant also reduces smearing of 
the Hubble flow by peculiar velocities (see, e.g., the numerical simulation in [17]). Therefore, 
the observations of the Hubble flows even at a few Mpc may reveal the presence of dark 
energy in the Universe [7, 10]. 

Now, we describe the differences between two our models. The main one consists in the 
spatial distribution of the inhomogeneities (gravitating sources). In the model with /C = 0, 
we were forced to introduce the radius Rq around a gravitating mass to ensure the finiteness 
of the gravitational potential at infinity. This radius is completely defined by the mass mo 
of the gravitating source and the cosmological parameters Hq and (at present time). 
Then, in such toy model, all inhomogeneities are surrounded by these spheres which we call 
the local gravity ones because outside of them the gravitational potential is identically equal 
to zero and we turn back to the unperturbed background Friedmann model with averaged 
ideal cosmological medium and the Hubble flows in it. In this model, the motion of a test 
body inside of the local gravity sphere occurs in full analogy with the Schwarzschild-de Sitter 
model where the acceleration is completely defined by the cosmological constant (without 
any admixture of dark matter); the zero-acceleration radius R* is at present time inside of 
the local gravity sphere: i?* < Rq if $7 a > and the dragging effect is absent for 

R> Rq. The spatial distribution of matter, where inhomogeneities are surrounded by empty 
spheres of local gravity and ideal cosmological medium is outside of these spheres, looks a 
little bit artificial. Therefore, we have considered the hyperbolic model QC = —1). Here, 
the inhomogeneities are distributed completely randomly. The finiteness of the gravitational 
potential for any number of gravitating sources follows naturally from the solution of the 
Poisson equation without any artificial cutoff. Outside of the zero-acceleration surface (with 
the radius Rh at present time), the dragging effect is not absent but goes asymptotically 
(rather fast) to zero. In contrast to the previous model, the gravitational interaction between 
distant inhomogeneities is not absent, although it is suppressed by cosmological expansion. 
We have also shown that the superposition principle for the nonrelativistic gravitational in- 
teraction works here. Therefore, our model gives a possibility to simulate the dynamical 
behavior of an arbitrary number of randomly distributed inhomogeneities inside of the cell 
of uniformity taking into account both the gravitational interaction between them and cos- 
mological expansion of the Universe. Some simple examples of such simulations are depicted 
on Figures 2-5. This is the main result of our paper. 
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